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ABSTRACT 


This report provides a brief description of a canonical set 
of equations which governs the behavior of an n-body spacecraft. 
General results are given for the case in which the spacecraft is 
modeled in terms of n rigid bodies connected by dissipative elas- 
tic joints. The final equations are free from constraint torques 
and involve only r variables (r is the number of degrees of free- 
dom of the system). An advantage which accompanies the elimina- 
tion of the constraint torques is a decrease in the computer run 
time (especially when n is large). 

Linearized models are obtained and are recast in the 
familiar form 

x(t) = F x(t) + G(x, t) u (t) 

where x is the state vector, u is the control vector, F is a con- 
stant matrix, and the matrix G depends on x and t. This form for 
the equations is particularly useful when modern control theory 
is used to arrive at a stochastic controller for a multi-hinged 
rigid-body spacecraft. 

The models provided in this report will be used in analyzing 
the cruise, the thrust vector control (TVC), and the articulation 
control (ARTC) modes associated with the Mariner Jupiter Mars 
(MJS T 77) spacecraft. Due to their generality, the models can 
also be conveniently used for analyzing a spacecraft appropriate 
for missions subsequent to the MJS ! 77 mission. 
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1 . 


INTRODUCTION 


This report provides a brief description of a particularly elegant 
formulation characterizing the rotational motion of a spacecraft idealized as 
a set of multi-hinged rigid bodies*, Assumptions made in this development 
include the following r 

(1) The spacecraft (S/C) can be adequately modeled as n hinged- rigid 
bodies connected by dissipative elastic joints, 

(Z) Chains of connected bodies do not form closed loops. 

(3) Only rotational motion is allowed at a joint. 

(4) There is a vector constraint torque orthogonal to the axis of 
rotation at a joint whenever the rotational motion has only one 
or two degrees of freedom. 

Considerable effort has been focused on the problem of obtaining the 
dynamical equations for an n-hinged rigid body spacecraft (e. g. , see Refs. 

[ 1] thru [5] ). The approach used in this analysis is based primarily on 
Refs. [1] and [ Z] . Likins (see Refs [4] and [5]) recently extended the method 
discussed in Ref. f Z] to appropriately account for flexible appendages. 

In this analysis, the appendages are considered rigid but viscoelastic 
joints are allowed. Effectively, the interaction torque existing at a joint 
connecting a purely rigid body and a flexible appendage is modeled as a spring- 
damper torque (that is, it is specified in terms of torsional stiffness and 
damping coefficients). This model is especially useful in the preliminary design 
and analysis of the spacecraft. 

The main objectives of this work are: 

(1) To obtain a linearized dynamical model of the n-hinged rigid-body 
spacecraft. 

(Z) To provide a brief description of some of the dynamical 
principles involved in the development of the S/C model. 
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The linearized dynamical model of the multi-hinged rigid-body 
spacecraft can be: 

(1) Used in a simulation to study the effects of interactions of the 
hinged members on thrust vector control (TVC) and articulation 
control (ARTC) performance, 

(2) Used in the development of a general stochastic controller 

(see Refs. [7] and [8] ) for a multi-hinged rigid-body spacecraft. 

1 . 1 Compendium of Report Contents 

In this section, a compendium of the section by section contents of the 
report are given. 

Section 2 gives the general nonlinear dynamical equations for an n-hinged 
rigid body spacecraft with r degrees-of-freedom. The main result of this 
section is given by the vector-matrix equation 

A co = L (1-1) 

where oo consists of the angular acceleration of the base body cLq and the relative 

angular acceleration components and the vector L consists of the forcing 

functions for the base body (Lq) and for the n - 1 remaining bodies (L 0 ). 

u R 

Section 3 provides a set of linearized dynamical equations for the 

n-hinged rigid-body spacecraft. An intermediate result of this section is 

given by the partitioned vector-matrix equation 


A 11 j A 12 


% 

A 21 1 A 22_ 


(I) 

L R J 


G 0 (y,t) 

G“cy,ir 


- 



(1-2) 


where the matrices A, B, K are constants, 0 is an appropriately dimensioned 
null matrix, and the matrices Gq and G^ depend on y (the relative angles) and 
t. The primary result of this section is the recasting of Eq. (1-2) in the form 


x(t) = F x(t) + G (x, t) u(t) 


(1-3) 
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where the state vector x consists of the attitude and angular velocity of the 
base body (0 and coq) and the relative angles (y) and angular velocities (co^) 
associated with the remaining n - 1 bodies. The control vector u(t) is related 
to the force Fq applied to the base body. 

Section 4 provides a linearized dynamical model for a specific topological 
arrangement of 5 -hinged rigid bodies. Expressions for the matrices G and F 
involved in the equation 

x = F x + G (x, t) u 
are obtained. 

Section 5 provides a linearized dynamical model for a 3-hinged rigid- 
body spacecraft. This result is a special case of that given in Section 4. 
Expressions for the linearized elements of the matrix A and the vector L 
appearing in the equation 

A oj = L 

are given. In addition, the expressions for the matrices F, G appearing in 

x = F x + G (x, t) u 
are given. 

Section 6 provides a simple linearized model for a single rigid body. A 
comparison of this result with the model for the n-hinged rigid-body spacecraft 
is made. 

Appendix A provides a discussion of the development of the dynamical 
equations for a multi-hinged rigid-body spacecraft (in terms of the constraint 
torques). The development as presented here is based on Ref. [1]. Inter- 
mediate results are tabulated to allow the reader to be cognizant of the origin 
of the various terms that are introduced,, 
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Appendix B provides a discussion of the development of a set of dynamical 
equations (free from constraint torques) for a multi-hinged rigid-body space- 
craft (see Ref. [2]). Intermediate results are tabulated for future reference. 

Appendix C provides a detailed description of an application of the 
general results of Appendix B. The explicit dynamical equations are obtained 
for a specific topological arrangement of 5-hinged rigid bodies. In Appendix C, 
the coordinate systems involved in the evaluation of the various terms are 
discussed. In addition, expressions for the elements of the matrix A and the 
vector L appearing in the equation 

A to = L 
are tabulated. 
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2 . 


DYNAMICAL EQUATIONS FOR AN n-HINGED RIGID-BODY 
SPACECRAFT 


In this section, the dynamical equations for an n-hinged rigid-body 
spacecraft are provided. Emphasis is placed on the procedure used to obtain 
the results rather than on a detailed and lengthy derivation of the results. 

Consider an n-hinged rigid-body spacecraft having r degrees of freedom. 
The number of scalar constraint torques for such a system is n^ = 3n - r. In 
Ref. [ 1] , Hooker and Margulies showed how to eliminate the constraint torque 
components so that 3n differential equations for the angular velocity components 
could be integrated (with the constraints being satisfied). In that technique, 
the calculation of the right hand sides (RHS) of these differential equations (DE) 
required solving a system of 3n + n^ = 6 n - r linear algebraic equations for 
the 3n angular velocity rates and the n^ constraint torques. 

In the present analysis, a set of r dynamical equations in which the 
constraint torques do not appear is given. The pivotal steps involved in obtaining 
this canonical set of equations involves (see Ref. [2]): 

(1) Recognizing that if the vector dynamical equations of all the bodies 
are summed, then the constraint torques cancel in pairs. 

(2) Noting that a vector constraint torque at a typical joint j can be 
isolated by summing the vector dynamical equations over all 
bodies that lie to one side of joint j (the constraint torques on 
this set of bodies all cancel in pairs, except for the one at joint j). 

(3) Observing that the constraint torque (isolated in step 2) at joint j 
is orthogonal to the gimbal axis at joint j. 

Effectively, 3 scalar equations result from the projection of the vector 
equations summed over all the bodies on to a suitable reference frame. More- 
over, r - 3 additional scalar equations result from the dot products of the r - 3 
gimbal axes and the constraint torques associated with these axes. The salient 
advantage associated with the elimination of the constraint torques is the 
accompanying reduction of the computer time required for integrating the 
equations (this is especially true for large n) ! 

In Ref. [3] , Fleischer describes a general computer simulation based on 
this technique. 
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Although more than a modicum of labor and more than a soupqon of effort 
were involved in the verification of the treatment given very succinctly by 
Hooker in Ref. [2], nevertheless, it is appropriate to keep the present dis- 
cussion brief. This brief treatment will allow attention to be focused on the 
underlying assumptions, on the interpretation of the results rather than on 
the detailed derivations of the multitude of lengthy equations ! 

The procedure used to arrive at the r scalar equations entails the 
following steps (see Table A - 1 ) : 

(1) Writing Newton 1 s and Euler* s equations for each body X. 

(2) Eliminating the unknown interaction force F. . 

j 

(3) Evaluating the term 



which represents the sum of the moments about the center of mass 
of body X. due to interaction forces F existing at joints j. 

j 

(4) Interpreting Euler's equations for body X (after using the results 

of step (3) ) as the equations for the augmented body X relative 

to its barycenter B x . 

K 

(5) Expressing the interaction moment acting at joint j on body X 

as a sum of a constraint torque hlf and a spring-damper torque 

, ,SD . K 

, i. e. , 

j 

,, ,,0 , SD 

M \. = M x. + M \. 

J J J 


The terms and symbols are defined in Table A-l as they are needed in the 
development. 
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(6) Recognizing that if the vector dynamical equations for the augmented 
bodies X are summed over all then the constraint torques cancel 
in pairs and consequently disappear, i. e. , 

EE M x = 0 

XeS jgJ^ J 

(7) Recognizing that the constraint torque at joint j acting on body \ can 
be isolated by summing over all bodies. to one side of joint j. 

(8) Recognizing that the gimbal axis g. is orthogonal to the constraint 

C J 

torque M. at joint j. 

Steps (1) thru (5) are discussed in turn in Appendix A; steps (6) thru (8) 
are discussed in Appendix B. 

2. 1 Compact Form for Dynamical Equations for an n-Hinged Rigid-Body 

Spacecraft 

In this section, the set of dynamical equations derived in Appendix B 
for an n-hinged rigid-body spacecraft are presented. In vector-matrix 
notation, the equations are: 



( 2 - 1 ) 


or 


A co — L 


In scalar form the equations are (see Table 2-1 for definitions of terms): 


r-3 

E 


a 00 * ^0 + 7 ■ a ov "Yv - k 


Ok Y k " ^0 - “X 

k=l \=0 


n-1 

E 


E, 


( 2 - 2 ) 
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r-3 


n-1 


a iO * w 


“o + 22 a ik \ = L i = % • 22 ‘i 


. E x , i = 1 , 2, . . . r-3 

K 


k= 1 


\ = 0 


whe re 


n-1 n-1 



\ = 0 |JL =0 


a 


Ok 



M- 


£ k[x % ®k 


a i0 





M- 


a 


ik 



e i\ % • 


E 


X. 


M X + D X x F X + 


E D, x [f + moo x (co x D )1 

XfJ- jx l M- M- M- |JA J 

(JL A 


- w x x \x 


+ 



r-3 


E k = E k 


22 % ■ 22 v ■ «k k 

[J. k=l 
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A comparison of Eqs. (2-1) and (2-2) reveals that; 

(1) Ajj is the 3x3 matrix representation of the operator • (where 

aQQ is a dyadic, and • represents the dot product). 

(2) A is the 3 x r - 3 matrix representation of the vectors a^* 

(3) ^22 r - 3 x r - 3 matrix representation of the scalars a^ 

(with i, k = 1, 2, ... r-3) 

Moreover, in Eqs. (2-1) and (2-2), 00 q represents the angular acceleration 
of the base body, represents the relative angular accelerations ( 7 ^) of the 
remaining n - 1 bodies, Lq is a 3 x 1 matrix; and is a r - 3 x 1 matrix. 
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Table 2-1. Definitions of Terms 


Item 

Definitions and Equations 

n 

The number of rigid bodies involved in spacecraft model. 

r 

The number of degrees - of-freedom of the system. 

d.. 

The vector from the barycenter of body X to the joint j 
of body X. 

D \| Jl 

D = D for all bodies p belonging to S . (the set of 

K\1 K j X J 

bodies 'connected to body X via joint j). 

D x 

Vector from barycenter of body X to c.m. of body X. 

□ x 

Inertia matrix relative to c. m. of body X. 

m x 

Mass of body X 

\x 

Augmented inertia matrix for body X relative to barycenter 
B X- r 

*xx = D x + m x ( D x' u - D x D x) + H m (1 ( D x^ u - % D x J 
L p A 


where U is unit dyadic 

^X(JL 

& = - m [d, • D , U - D X D, ] 

Xp L Xp pX pX Xp J 

§k 

Gimbal axes, k = 1, 2, ... r - 3; g^_ is a unit vector 

€ kfj. 

= 1 if gimbal axis g^_ is between body p and body o, 
otherwise € = 0; € specifies bodies p which sense the 

rotation Y k g k * 

F x 

Vector representing external force applied to body X. 

M x 

Vector representing external moment applied to body X. 

c x 

H- 

Direction cosine matrix transforming coordinates of body p 
to coordinates of body X. 

A ,SD 
T» . = M. 

Xj Xj 

Vector representing spring-damper interaction torque on 
body X at joint j. 

J x 

Set of labels for joints j belonging to body X. 
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3 . 


LINEARIZED SET OF r DYNAMICAL EQUATIONS FOR AN n-HINGED 
RIGID-BODY SPACECRAFT 


In this section, a linearized set of r dynamical equations for an n-hinged 

rigid-body spacecraft is provided. Linearization is accomplished by retaining 

only terms of first order in coq, y^ an< ^ their derivatives in the solution (i, e. , 

products of gOq and y^ with k = l, 2 , ... r -3 and their derivatives are neglected). 

In addition, it is as sumed that y, (with k = 1 , 2, ...r-3) and 0 . (with i = 1 , 2, 3 ) 

ic 1 

are small angles -- hence the direction cosine matrices take a particularly 
simple form. 


3 . 1 Direction Cosine Matrices 

Typical direction cosine matrices for the linearized case become (see, 
e. g. , Fig. 4-1) 

^2 ^ 

C 1 = E - v 2 b 2 

C t = E - V 4 S 4 < 3 -') 

C 0 = E - V, 8j 

C 0 - C 3 Cq = [E - y 2 B 2 ][E - Yl g,J*E- V, g, - V 2 g 2 
c 0 = 1 E - Y3 g 3 ] 

C 0 = c 3 c 0 = [E - V 4 g 4 ] [E - y 3 g 3 ] * E - y 3 g 3 - Y 4 g 4 
0 

C AT = E - 0 = E - 
N 


o 

1 

CD 

1 

6 Z 


1 

iH 

CD 

rvl 

CD 

i 

1— 1 
CD 

0 

- 6 3 

= 

1 

CD 

i — 1 

l 

6 3 

I 

CD 

0 3 

o 


6 2 

- 0 3 

1 


Products such as u>qY and yy are neglected in the linearization in this report; 
such terms can be retained and included in the forcing function L if it is 
desirable ! 
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where the vector 0 consists of ordered rotations 0^, 0^, 0^, E is a 3 x 3 
identity matrix, and - over a vector represents the matrix representation of 
the cross-product operation. 


Relationship Between Attitude and Angular Velocity of Base Bod^ 


The relationship between the attitude rate and the angular velocity of 
the base body 0 becomes (for the linearized case) 



when the small angle assumption is used and in addition products of 0^ and 
are neglected. 


3.3 Evaluation of Elements for Linear Case 

In this section the elements a^ m are evaluated for the linear case. 
Recognizing that products of y^. and y^ and y^ and d>g can be frequently neglected 
(for k = 1, 2, ... r-3), it is clear that only those portions of that are not 

functions of y^ are to be retained. Recall from Eqs. (2-1) and (2-2) that the 
a j0m. ,S are multipliers of cLq and Effectively, this implies that the 

direction cosine matrices C^_ (withk = l, 2, . . . r-3) appearing in the expressions 
for the a ^ m * s can approximated by identity matrices. The matrix A which 
is composed of the elements a ^ m then becomes a constant. 

3 . 4 Evaluation of Forcing Function L for Linearized Case 

In this section, the terms involved in the evaluation of the forcing 
function L are provided for the linearized case. Recall that L is defined 
according to (see Eq. (2-1)) 


Ad) = L 


(3-3) 


where d) consists of d)g and ob-^ and L consists of Lg and L^ t 


12 


JPL Technical Memorandum 33-613 



Recall, too, that L and L , the components of L, are given by 

0 K 

L 0 - E E x 

X 


0-4) 

81 ’ ? ' 1X ^ I 


e Z • E *2X K 


| g r-3 * € r-3,\ E X ^ 

The linearized versions of E and E reduce to 

K K 


E x = M X + D X X F X + 


E 


D, X C X 

Xq. (JL 


F + 

V- 



T 


>0 


E 



E. 


(3-5) 


and consequently, the linearized versions of Lq and reduce to 


l o= E' 

M x + D x x F x + 

/ D x F + 

Z— J Xq (-L fj. 

E I 

(3-6) 

X 


q. A 

j £j x ) 
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n-1 


[ 1 • X ‘>K D 


xo x c o F o + g l • / Jl\ 


\=1 


\ j 6 J 


S'N 


R 


n-1 


Si ’ / £ D n x C? F. + g 0 

& Z / i 2X XO 0 0 6 2 

X = 1 


X *2X X T 


\ j€j. 


X . 
J 


X 


n-1 


g. ^ ‘ ^ «« i i D xn x C n F n + S 


r-3 ^ C r-3,X XO A 0 X 0 s r-3 "r-3,X ' Xj 

X= 1 X jeJ 


XI e r-3,X X 


3. 5 Compact Form the for Linearized Set of r Dynamical Equations 

In this section, the linearized set of r dynamical equations for an 
n-hinged rigid-body spacecraft are expressed in compact form. 

First, the term 

vEvEv <3 - 8) 

x ,., x ^ 

is examined. As pointed out in Appendix B, summing the dynamical equations 
over bodies X., which are connected beyond gimbal axis g^ relative to the base 
body, isolates the interaction moment at joint i on body \. This implies that 

e iX X T x. 5 g i • T ii = g i • 

x je j x J 


1 1 


By' 

li 


g i 


(3-9) 
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whe re 


11 


K.v.g. - B.v.g. 

i Y i & i i Y i s i 


In Eq. (3-9), and are the stiffness and damping coefficients associated 
with joint i. Substitution of Eq, (3-9) into Eq. (3-?) yields 


L A = 


n- 1 

M o + D o x F o + Xj D X0 x c o f o 

\=1 


(3-10) 


I 


L 


R 


«i • E 'i\ D xo xC o F o - < K 1^1 + Vl 


g 2 ’ 


E' 2 X° 


X0 x C 0 F o - < K 2 v 2 + B 2 v 2 > 


’r-3 


E 


e r-3,\ D \0 X C 0 F 0 


(K ~ y o + B ^y o ) 
r-3 * r-3 r-3 1 r-3 


Equation (3-10) can be written as 


^0 + F 0 + C 


F q = G 0 ( V ,t) Fq 


(3-11) 


L R = G R^ jt) F q - K V - B Y = G R (Y,t) F 0 - K v - B 7 


16 


JPL Technical Memorandum 33-613 



where 


Gq(y, 1;) is a 3 x 3 matrix 
G-n>(Y,t) is a r-3 x 3 matrix 

x\ 

K, B are r - 3 x r - 3 diagonal matrices 

~ over a term represents the skew symmetric matrix representation 
of the cross product (e. g. , x = ^-\q) 

Y, y are r - 3 x 1 matrices consisting of elements Y^> 
for k=l, 2, . . . r-3 

Equation (3-11) can thus be written as 



or 


L = G F 0 



(3-12) 


where 0 is an appropriately dimensioned null matrix. 

Collecting the results of Sections 3-3, 3-4 and the results of Eq. (3-12), 
it follows that a set of r linearized dynamical equations for an n-hinged rigid- 
body spacecraft is given by 


hl_ 

;_ A ji 

/ W ° \ 

Go 

TP 4- 

0 

■v / x 

0 

A 21 

1 A 22. 

\ W R / 

. . 

F o + 

-K 

Y + 

-B 


or 


A to = G(y,t) F q + 



where A, K, B are constant matrices and G depends on y and t. 


(3-13) 
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3 . 6 Linear Dynamical Model for n-Hinged Rigid-Body 
Spacecraft (State Equations) 

In this section, the linear model developed in Section 3. 5 is cast in a form 
suitable for use in modern control theory. Essentially, the state equations 
are sought. As seen in Section 3. 2, 



and 


Y = co R 

where the vectors 0 and coq are the attitude and angular velocity of the base 
body relative to an inertial frame and the vectors y and oo^ represent the relative 
attitude and angular velocities y^> y^ for k = 1, 2, ... r-3. The state can 
thus be defined as the 2r x 1 vector 


x 


\ 

w o 1 

i 

V “r / 


The differential equations for 0 and y are given above and those for 
cOq and can be obtained from Eq. (3-13). 

Manipulation of Eq. (3-13) yields 


A 11“0 + A 12“R = G o'^ t)F o (3 - 14) 

A 21 "0 + A 22 “R = G R ( ''> t) F 0 - K V - B-i 
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As discussed in Appendix C, Eq. (3-14) can be written as 


W 0 


A 1 1 ‘ A 12 A 22 A 21 


-1 


F 0 


- Aj l - A 1Z A 22 A 21 A 12 A 22 ( G R^’ ^ F 0 " ) 


-1 


w 


j R = A 22 I G R (V»t) F n - Ky - B\ - A ?1 (I) 


21 


Redefining the bracketed matrix as a, it follows that 

“0 = a ' 1 F 0 - a - 1 A 12 A-* {g r (y> t) F 0 - K V - By } 

W R = A 22 | G R^’^ F 0 - K ^ - A 21 w 0 

In vector-matrix rotation, the state equations become 


(3-15) 


(3-16) 
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where E is a 3 x 3 identity matrix and 0 is an appropriately dimensioned 
null matrix. Note that Eq. (3-17) has the same form as 

x = F x + G(x, t) u (3-18) 

where x is the state and u is the control variables. It is important to note that 
F is a constant and G depends on x and t! Equation (3-17) represents the 
primary result of this report. 

It is immediately recognized that the solution to Eq. (3-18) can be 
written as 


x(t) = <(>(t, 0) x(0) + 


<f>(t) 4> _1 (t) G(x, t) u(t) dT 


(3-19) 


where x(0) is the initial state and cf>(t, t) is the transition matrix. The transition 
matrix can be computed from the matrix differential equation 


i = F cj> 


(3-20) 


with 


<M0) = I 


or it can be computed analytically if n is small. 
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4. LINEARIZED DYNAMICAL MODEL FOR A 5 -HINGED 
RIGID-BODY SPACECRAFT 

In this section, the general linearized dynamical model developed in 
Section 3 is used to obtain expressions for the elements of the constant matrix A 
and the vector L for the case in which n = 5 (see Fig. 4-l)„ The expressions 
for the elements of the matrix A are summarized in Table 4-1 and those for the 
elements of the vector L are summarized in Table 4-2. Note that in Table 4-1 
the subscripts appearing to the right of a term in parentheses refer to the coor- 
dinate systems in which the terms are computed or expressed. 


The 

7x7 matrix A is 

given by 



\ a 00 a 01 a 02 
\ 

a 03 

a 04 

A = 

N ^ll a 12 

\ 

a 13 

a 14 


\22 

\ 

a 23 

a 24 


symmetric \ a 33 a 34 

X \ 

\ 44 


The partitioned matrix A (involving the 3x3 matrix A^, the 3x4 matrix A^, 
the 4x3 matrix A^, and the 4x4 matrix A^) is given by 





nn 

° 

1 o 

1 

_L 

a oi 

a 02 

a 03 

a 04 




a io 

I 

l 

i 

a n 

a l2 

a 13 

a 14 

fnj 

L A _>i 


a 20 

l 

1 

i 

a 21 

a 22 

a 23 

a 24 

A 21 | 

A 22 


a 3 0 

1 

1 

i 

a 31 

a 32 

a 3 3 

a 34 




a 40 

l 

1 

a 41 

a 42 

a 43 

a 44 
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The 3x3 time- varying matrix Gq is given by 


G 0 = ? 0 +f> 0 + ( C ? D 11 +C 2 D 22 + C ? D 33 + C ! D 44)'' 

The 4x3 time -varying matrix G 0 is given by 

XV 

^( D ll +C 2 D 2 z)’ c 0 

6 22 C 0 

*3 (d 33 4 D 44 [ 

T fS r^ 4 
g . D. . C n 

6 4 44 0 


(4-3) 


(4-4) 


The constant diagonal matrices K and B are given by 


K, 


B, 

1 


1 

K. 


B_ 

2 


z 


; b = 


K_ 


B 0 

3 


3 

K . 


B , 

4 


4 


(4-5) 


and the constant 3x3 matrix a is given by 


A 11 " A 12 A 22 A 21 


(4-6) 
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Hence, all matrices involved in the generic form [see Eqs. (3-17) and (3-18 
x = Fx + G (x, t) u 

have been defined! Note that the symbol - appearing as 

( f 

indicates that the tilde is to be applied to the resultant expression within the 
parentheses. 



Figure 4-1. Pictorial Sketch of a 5-Hinged Rigid-Body Spacecraft 
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Table 4-1. Expressions for Linearized Elements a . a 

(i, m = 0, 1, ... 4) for a 5 -Hinged ~' in 

Rigid-Body Spacecraft 


Equations 


[*00 + *01 + *02 + *03 + *04 ^ 0 
+ [* 10 + ®n + *12 + *13 + *14 ] 1 

+ f*20 + *21 + *22 + *23 + * 24^2 

+ ^*30 + *31 + *32 + *33 + * 34 ] 3 

+ [*40 + *41 + *42 + *43 + * 44] 4 


| [ * Ql + *02 ^ 0 + ^ *1 1 + *12 b + ^*21 + *22 ^ 2 

+ [®31 + *32 V (*41 + *42 ^ 4 ] ' < g l>l 


^* 02^0 + ^* 12^1 + ^* 22^2 + f ® 32^3 + ^* 42^4 ' ( g 2*2 


1*03 + *04 V [®13 + *14V [ ®23 + ®24 ] 2 


+ [* 33 + *34)3 + [*43 + *44)4 


t®04l 0 + [ ®14ll + t®24i 2 + t*34^3 + '* 44)4 ' < g 4>4 


(g l>l ' t*n + ®12ll + 1*21 + *22^2 ' < g l>l 


(g^j * [* 12 ^ + [* 22^2 * [ g 2^2 


[ g i^l ' [ *13 + *14^1 + ^ *23 + *24-^2 * * g 3^3 
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Table 4-1. Expressions for Linearized Elements a 
(1!, m = 0,1, ...4) for a 5 -Hinged 
Rigid-Body Spacecraft (contd) 


Item 

Equations 

a 14 

( g j)i 

j ^14^ 1 + ^24] 2 

• (g 4 ) 4 

a 22 

( g 2^2 * ^ 22^2 * ^ g 2^2 

a 23 

(g 2 } 2 ' t*23 + $ 24^2 ' (g 3>3 

a 24 

(g 2 ) 2 ' [*24^2 ' (g 4 } 4 

a 33 

(g 3 ) 3 • 

[$33 +^ 34 ] 3 + [$ 43 + $ 44 ] 4 * < 83)3 

a 34 

( g 3 ) 3 

[$ 34 ] 3 + t $ 44 ] 4 

' (g 4 ) 4 

a 44 

(g 4 ) 4 • [ $ 44 ] 4 ‘ (g 4 ) 4 

^0 

(a Qk ) T , k = 1 , 2, ... 4 

a £m 

^ a m£ ^ T ’ k = 2 > * * • 4> 1 ^ k 
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Table 4-2. Linearized Forcing Function L for a 5-Hinged 
Rigid-Body Spacecraft (External Forces and Moments 
Applied Only to Base Body) 


Item 

Equations 

L o 

M 0 + D 0 x F 0 

+ ( C ? D 11 tC 2 D 2Z + C 3 D 33 + C 4 D 44) xF 0 

L 1 

• ( D ll tC 2 D 2 2 ) xC 0 F 0- K 1 V! - B; Y, 

L 2 

*2- ( D 22 xC 0 F o) - K 2 V2- B 2^2 

L 3 

S 3 ' ( °33 +C 4 D 4 4 ) xC 0 F 0 - K 3 V 3 - B 3 V 3 

L 4 

g 4 ' ( D 44 xC 0 F o) ' K 4 V 4 - B 4' } 4 
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5 . 


LINEARIZED DYNAMICAL MODEL FOR A 3 -HINGED RIGID -BODY 
SPACECRAFT 


In this section, the linearized dynamical model for a 3 -hinged rigid-body 
spacecraft is given/ This result is a special case of that given in Section 4 and 
is obtained by defining the new base body to include bodies 0, 3, 4 of Fig. 4-1. 
Figure 5-1 shows a pictorial sketch of the resulting 3 -body system. 

The expressions for the elements of the 5x5 constant matrix A and the 
vector L are given in Tables 5-1 and 5-2, respectively. 

The 5x5 constant matrix A is given by 


00 

a 0l 

a 02 

10 

a ll 

a 12 

20 

a 21 

a 22 


The partitioned matrix A (involving the 3x3 matrix A^, the 3x2 matrix A^, 
the 2x3 matrix A^ and the 2x2 matrix A^) is given by 


A 

A 


21 I 



P 1 
o 
o 

p 

O 

p- * 

a 02 

^ 1 

o 

1 

a n 

I fM 

i • 

1 

1 

a 20 1 

a 21 

a 22 


The 3x3 time -varying matrix Gq is given by 


G o - i 0 +D 0 + 


( C ? D 11 +C ° D 2z) 


(5-2) 


(5-3) 


t 


This model was used to obtain the results given in Ref (6). 
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The 2x3 time- varying matrix (Gp is given by 


R 


*T ( D 11 +C 2 D 22 ) 


T fC r 2 

L g 2 D 22 C 0 


(5-4) 


The constant 2x2 diagonal matrices K and B are given by 


K = 


K, 



B 


B. 



(5-5) 


and the constant 3 x 3 matrix a is given by 


a - A n A j2 A 22 A 21 J 


(5-6) 


Hence, all the terms needed in the generic form of the state equations 
given in Eq. (3-17) and (3-18) have been specified. 
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Table 5 - 1 . Expressions for Linearized Elements aj0 m (f,m = 0 , 1 , 2 ) for 

a 3 -Hinged Rigid- Body Spacecraft 


Item 

Equations 

o 

o 

[®00 + *01 + *02) 0 + f *10 + *1 1 + *12 'l 
+ [* 2 „ + * 21 + * 22 ] 2 

a 01 

{[*01 + *02)0 + [®11 + *12^1 + [®21 + *22^2 J* <*1»1 

a 02 

| ^ ^0 2 0 + f ®12b + ' * ^ g 2^2 

a ll 

B 

[$1 1 + $1 2 ] 1 + f$ 21 + $ 22 J 2 * g l*l 

a 12 

B 

( [*12]l + t *22^2} • <«2>2 

a 22 

^ g 2^2 ' 


f $ 22^ 2 ' * * g 2^2 

a io 

r 

( a oi )T 

a 20 

( a o 2 ) T 

a 21 

< a 12) T 
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Table 5-2. Linearized Forcing Function L for a 3-Hinged Rigid-Body 
Spacecraft (External Forces and Moments Applied Only to Base Body) 















6, 


LINEARIZED SINGLE RIGID-BODY SPACECRAFT MODEL 


In this section, the linearized model for a single rigid-body model of a 
S/C is given. In this case, the generalized results for a n-hinged rigid body 
are not directly applicable. Instead, the linearized version of Euler's equa- 
tions for a rigid body are used. That is, the state equations are 


il 

•^1 

' - 0 J _ E - 1 

o 

•3 

i 

o 

o 



+ 



where 

E is a 3 x 3 identity matrix 

Mq is the external moment applied to the S/C 
I is the 3x3 inertia matrix of the S/C 
0 is a 3 x 3 null matrix 


( 6 - 1 ) 


It is interesting to note that the linearized form for an n-hinged rigid-body 
spacecraft can still be used for a single rigid body model of the S/C, This is 
achieved by eliminating the vectors y and co^ from the state and appropriately 
interpreting the results. Comparing Eqs. (3-17) and (6-1), it follows that for 
the single rigid body model 


a = I (the inertia matrix) 


G 


0 


t 


0 


and y, , G^ , K, B do not appear (they are deleted from the general r'esult). 

Xx Xv 
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APPENDIX A 


DYNAMICAL EQUATIONS (IN TERMS OF CONSTRAINT TORQUES) 

FOR A SPACECRAFT IDEALIZED AS n-HINGED RIGID BODIES 

The primary result of this appendix is a compact version of the dynamical 
equations for a spacecraft idealized as n-hinged rigid bodies; the equations are 
written in terms of the constraint torques. The pivotal steps involved in the 

,A. 

derivation are briefly described and the intermediate results are recorded 
(see Table A-l). 

The steps involved in obtaining the desired result include: 

(1) Writing Newtons and Euler’s equations for each body X. 

(Z) Eliminating the unknown interaction force F 

X . 

J 

(3) Evaluating the term 



(The sum of the moments about the center of mass of body X due 
to interaction forces F\. existing at joints j). 

(4) Interpreting the results as the equations describing the motion of 
the augmented body X relative to its barycenter. 

A-l Newton’s and Euler’s Equations 

The development of the equations of motion for an n-hinged rigid-body 
spacecraft begins with Newton’s and Euler’s equations written for each body X. 
That is, for all Xe S: 


'The intermediate results are important since they indicate the origin of 
various terms that appear in the final result. 
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F X + ^ F X. ■ m X (fe + Px' = “X F c. 

r l X 

J <J X 

(A - 1 ) 

Ox“x + "x^x • “x = M x + 'll M x. + Y! C X, x F x. 

'•x j A 3 3 


where 

0 is the inertia dyadic for body V relative to its center of mass 

m. is the mass of body X 

F. is the external force applied to body X 
X 

oo. , w are the angular velocity and angular acceleration of body X 
X X 

relative to an inertial frame 
is the interaction force acting on body X at joint j 

j 

R is the position vector from the origin of an inertial reference frame to 
the point "O" ("O" is the origin of the S/C coordinate system in 
the undeformed state) 

P x is the position vector from the point "0" to the center of mass of 
body X 

R is the position vector from the origin of the inertial reference 
C X 

frame to the center of mass of body X 

R c is the position vector from the origin of the inertial reference frame 
to the composite mass center of the system 

is the external moment applied to body X relative to its center 
of mas s 

is the interaction moment acting on body X at joint j 

j 

= d^ is the position vector from the center of mass of body X 

j j 

to joint j 
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S is the set of all labels for bodies \ ^that is, S = jo, 1, . . ■ n - l} j 
is the set of labels for all joints associated with body \ 

A- 2 Elimination of the Interaction Forces 

In order to obtain a solution forco^ from Eq. (A-l) the interaction forces 
F\. are first eliminated. This is accomplished by recognizing that 



(A- 2) 


where S\ . 


is the set of those bodies connected to body \ at joint j. 


That is, by summing the interaction forces acting on body \ over all 
bodies p belonging to the chain of bodies connected to X. at joint j, the inter- 
action force can be isolated. Using Eq. (A-l), it follows that 




m 


R ) 

c 


(A-3) 


The form of F^. given by Eq. (A-3) is next substituted into Eq. (A-l) to 
eliminate F^. from the equations, that is 


D . • w-. + w. x 0. • 


x “x + “x x U X ' “X = M x + 


E m x, + E c x, x E <w*c > 


c x. x 

j tJ X j J eJ x j 


A-3 Evaluation of the term 


E 


c x. x F x, 

J 3 


j £j x 

In this section, attention is focused on the term 


(J.€S. 


X. 

J 


!*• 

(A-4) 


E 


c x. x F x. : 

j eJ v J J j«j 


E C. xV (F - m R ) 

X [J. |T c 
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It can be easily shown that 


E C, xF, = D. x F. + / D. x F - / m C x (R - R ) 

X. X. XX / j K|jl (jl (J. Xfi c c' 

: , t J J , , 4\ 1 1 -f\ l 1 


where 


D, 


E 

\x A 


m 


H* 


m 




D. 


\\x 


D, + C 


Xp 


If the point masses (the mass of all bodies attached to X. via joint j) located 

at joints j are augmented to the mass of body X, the augmented body X results. 
The baryc enter is defined as the new center of mass of the augmented body 
X, Physically, D is the vector from the barycenter B. of the augmented body 
X to the original center of mass of body X, is the vector from B. to joint j. 

j *■ 

Note that 



D v for all [jl tS. 



Next the term 


-E 

(-*• A 


m 

H- 



x (R 

c 

M- 



is examined. Using the facts that 


R - R = D + 
ecu. 

H- ^ 


E 


D 


V[x 
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0 


D = D - C 

|X (JLV [XV 


m. D. + / m D 

X X. / j (x Xp 

fx A 


it can be shown that 


L m C. x (R - R ) = - X • w - w x X . x> 

|X X(x c c XXX XX 

[x A |1 

+ m / D x D . 

/ j Xp )X X 

ix A 

where 

X = m. D . 2 U - ni h D. + V m ( D . 2 U - D D ) 
XX X X X / j (J, \ X[i \[X \\i ) 

H A 

U is an identity dyadic 

and are dyadic s. 


(A-6) 


Note that physically represents the inertia dyadic which must be added to 
[]^ to yield the inertia dyadic of the augmented body X relative to its barycenter 


B 


X* 


The term 


ix A 


\ x D x 
X|x (X A. 


can be expanded to yield 


JPL Technical Memorandum 33-613 


37 



m / D. xD , = m ) D. x xD +co x(oj x D ) 1 (A- 7) 

/ j Xp pX / J Xp 1_ |JL pX p p |JL X J 

|JL ^X. (J- A 

since D is assumed fixed in body p (the rigid body assumption). Collecting 

[JL A. 

the results of Eqs. (A-5) thru (A-7), it follows that 

E C. x F x = D N x F x + / D x F - X. • -w. x X. • w. 

X. A.. X. X / J Xp p X X X X X 

jeJ K J J pA 


+ m 


E 

p A 


D 


Xp 


x ^co,, x D ^ + oj _ x (u>" x D 


P 


pX 


P 


p X 


>] 


(A-8) 


If Eq. (A-8) is substituted in Eq. (A-l), an elegant result is obtained -- the 
equations governing the behavior of the augmented body X realtive to its bary- 
center B x ! 


A -4 Compact Form for the Rotational Equations Characterizing Body X 
Written Relative to its Barycenter 

Use of the results of Eqs (A-l) and (A-8) allows an interesting interpre- 
tation of the equations of motion for body X to be made. Substitition of Eq. (A-8) 
into Eq. (A-l) yields 



+ D x F + / D x F - X • & 

X X / j Xp p X X 

p A 


co N x X x • oa + m 


"X ~ “X ' ~x ' “"Xu °Xp x (cL p x D p^ 
p A 


+ m / D xLo x (u> xD.) 

Z — i >^p L p h- p^ J 
p A 


(A- 9) 
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Manipulation of Eq. (A-9) yields 


x x x • w v = M v + D x x F, + ^ D, x F 


"XX X X "XX X 


X X X Xp “ p 

P A 


+ m XI V x [ i 

|iA 


o>. . xD N + 0 J,, x( Wii xD ) 


J-/ i uj w v (jj 

p pX P P 


X 


M 


X. 

j« j K ■> 


(A- 1 0) 


where is the inertia dyadic of the augmented body relative to the baryc enter 

B X- 

Equation (A- 10) is the main result given in Ref j^lj . Briefly, Eq, (A- 10) 
implies that the rotational motion of body X of an n-body system can be 
obtained by 

(1) First forming the augmented body X by adjoining the masses m^ 

occurring at joints j belonging to J to the mass of body X ( m * ).^ 

X X 

(2) Determining the inertia dyadic of the augmented body X according 
to 




= o 


*XX U X + m x D x u " “x U X U X ' /__J "*p L“Xp “ "Xp "Xp 

P A 


D. D. + 


X % K u - ■w. 


(3) 


Considering the terms M. + D x F + / D 

XXX 

P A 


. x F as the 
Xp p 


external moment applied to the augmented body X relative to its 
barycenter . 

(4) Considering the terms 


m 


X V X % X V’ + m X V X K X % X V>] 


P A p A 

to be due to "inertial" forces. 
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Consequently, a particularly elegant and appealing result is obtained! In fact, 
the form of Eq. (A- 10) is reminiscent of that appropriate for a purely rigid 
body relative to an arbitrary point P, viz. , 


• ob + oo x • oo = 


m d x R.^ 
pc P 


(A - 1 1 ) 


where 


is the inertia dyadic relative to the point P 



is the vector from the inertial frame to the point P 

is the vector originating at P and terminating at the center of 
mass c 


oo, ob are the angular velocity and angular acceleration relative to an 
inertial frame. 


A compact form for Eq. (A- 10) can be obtained by manipulating its terms. 
First, the term 


m 


H % x (l V xI Vx> 

[j. A 


is examined. It is clear that 


D x (w x D ) = - D x (D x w ) 
X.|J. |Ji |a X. X.[x |jA |j. 


(A- 12) 


and that 


V x 


(D xcb ) = - f D . • D U-D.D, I • w 
pX p L pX Xp pX Xp J fji 


where U is the unit dyadic and is a dyadic (note that the vectors 

and are not fixed in the same bodies). 
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Substitution of Eq. (A-1Z) into Eq. (A-10) yields 


\\ ■ + % 


M* 



or 


E 

| 0 , e S 


$ '03 = E + 

XfJL fJL X 


E 

jeJx 


M. 


X. 

J 


(A- 1 3 ) 


where 


V 

\x 


- m 


[V D \|x U D ^x\ D X,x] ' ^ K 


Ox + “V D x • D x u ' D x D x> + E W ' V u - W 


|a A 


E, = M. 


+ D X x F X + E V X F |X + E % X [ m °V X '“n X V 1 ] 

fj. i X H A 


- “x x \x ’ w x + 7 M - 


E 


jej, 


_SD 

X. 

J 


A/r Ayr C j a/T^D 

M x. = M x. + M x. 


J J 


J 


In Eq. (A- 13), note that the interaction moment is assumed to 

consist of a constraint moment and a spring-damper moment. Recall that 
modeling the interaction moment in terms of a spring-damper is a consequence 
of the assumption that the joint is dissipative and elastic. 
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Table A-l. Rotational Equations for an n-Hinged Rigid-Body Spacecraft (In Terms 

of Constraint Torques) 


Equations 


Newton' s and 
Euler's equations 
for body \ 


F \. = m \ (i * + "v = m x v 




h S M X. + 5Z c x. xF x 


m is mass of body X 

c^ is vector from c. m. of body X to joint j; 

is set of labels for joints on body X; 

S\j is set of these bodies connected to body X at joint j; 
S is set of labels of the ensemble of bodies; 

is set of labels for joints on body X; 

F^ is external force on body X; 

F\ is the interaction force on body X due to joint j; 

j 

is external torque on body X 

is interaction torque on body X due to joint j; 
is the inertia dyadic for body X 


Elimination of 
unknown inter- 
action forces 


(F - m R ) 


Summing the interaction forces F^ over all bodies 
belonging to isolates F^ 

j j 


Evaluation of 


C X. xF X. 
J J 


c x. x F x. 


D X xF \ + . 


E m c. x (R - R ) 

p Xp C c 


R c is vector from origin of Newtonian frame to composite 
mass center; 

barycenter of body X is the new center of mass 
obtained by augmenting the point masses 

m^ (the mass of all bodies attached to X via joint j) 
located at joints j to the mass of body X; 


-1 

m m c x 
p Xp 


is the vector from barycenter to c. m. (c^) of body X 
is vector from to joint j; 


= D for all p* S. 






[il 


Table A-l. Rotational Equations for an n-Hinged Rigid-Body Spacecraft (In Terms) 

of Constraint Torques) (contd) 
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Resulting 
rotational equa- 
tions for body \ 
after evaluating 


c \. x F x. 

J J 


U ■ u x + w x x U * ui x + 7 M. + D. x F, 


"X X \ X X 


X. X X 


E D x F - X -u -U xX, 'a 
Xp p X X X X X 


Obtained by substituting expression for la 


into Euler’s equation for body X; the term M + D x F. 
is the moment relative to the barycenter due to 
external force F v ; the term D* + X is the inertia 
dyadic of the augmented body X relative to the bary- 


17 D x xTcj xD 
Z— < L M- I 

P A 


x(w xD.) 

M- ^ h 1 ' p P>^ J 


{ U x + x x } ‘ “x +w x x ( ll x + x x> * w x 


+ D x x F v + 


E vE 

j € J x J MX 


x F + m 


(u x D 

\ V- ^ 


x x (oj x D x ) 

X p p |JL X 


Compact form 


XI v ■ “ 

P € S 




E x = M. + D x 
K \ K 


M\. is considered the sum of an unknown constraint torque 
= ^ and a spring-damper torque (M^ = M\ + ); note that 

jeJ^ ^ the vectors used in defining ^ are not fixecl in the 

same bodies; the following relationship involving a cross 
P + \ ' j*) x Tf + m x( xD )] product and the dot product of a dyadic and a vector is 

X Z—rf Xp‘ L M- p '°V X pX J j , , . . 

^ used in obtaining 


x \\ • E 


3 


V = - m [V’ 


D. U - D x D x ] for p ^ X 
Xp pX Xp_ 


E + X x = u x + 


m X u -W + £ 


m (D“ U 
p v Xp 


Ax(BxC) = A.CB-A'BC = -[A'BU-BA] • c 


D. D ) 
Xp Xp 







APPENDIX B 


A SET OF r DYNAMICAL EQUATIONS (FREE FROM CONSTRAINT 
TORQUES) FOR A SPACECRAFT IDEALIZED AS 
n-HINGED RIGID BODIES 

In this Appendix, a set of r dynamical equations free from constraint 
torques is given (see Table B-l). As noted in Section 1, the procedure for 
eliminating the constraint torque involves: 

(1) Recognizing that the sum of the vector dynamical equations for 
each X over all bodies contains no constraint torques (they cancel 
in pairs according to Newton's Third Law). 

(2) Recognizing that a vector constraint torque at a typical joint j 

can be isolated by summing the vector dynamical equations over all 
bodies that lie to one side of joint j (the constraint torques on this 
set of bodies all cancel in pairs, except for the one at joint j). 

(3) Recognizing that the constraint torque at joint j is normal to the 
gimbal axis (axis of rotation) at joint j (this follows from the 
definition of a constraint). 

B-l Summing the Dynamical Equations for Body X Over All Bodies 

In this section, the equation obtained by summing the dynamical equations 
for body X over all the bodies of the system is examined. This result is given 
by (using Eq. (A-13)) 


EE vvE e ^EE m : < b -‘> 

Xfji X XjeJ^J 

As noted previously, however, the constraint torques vanish when they are 
summed over all j € and all Xe S, i. e . , 

E E m x. = ° (b - 2 > 

x 3 
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Furthermore, if there is a chain of bodies connecting a particular body p to 
the so-called base body designated by 0, then 


00 = G0 n + 

p 0 


r-3 

E« 

k= 1 


kfi. Y k g k 


(B-3; 


where 


oo q is the angular velocity of the base body 

g is the unit vector representing the gimbal axis at joint k (there 
are r-3 gimbal axes in all) 

y v is the relative angular rotation of the two bodies connected at joint k 

k g k 


is 1 if body p senses the relative rate an< ^ ® if body p does 


not sense it. 


It is clear that if joint k lies between body 0 and body p then = 1, 

Kli 

otherwise = 0. Using the result of Eq. (B-3) it follows that 


r-3 


r-3 


°V = "o + 2-J 'kp. V k 6 k + z_, 'k^ V k 8 k 

k= 1 k=l 


V n , g v + € 


(B-4) 


Substitution of Eq. (B-4) into Eq. (B-l), yields 

r-3 

EE 


k P 


■E *>-E 

X. u k= 1 



J. - -f 

E 

k=l 


€ k(J. g k Y k 


€ k|j. Y k g k 


(B-5) 
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^ 1 *kp ^k[JL 

k=l \ p 





In particular, the operator (dyadic) 


E Ev 

X p 


€ 


kp 


(B-8) 


occurring in the RHS of Eq. (B-8) is the same as that used in defining the 
term a^^ as expected. 

B - 2 Obtaining r-3 Additional Scalar Equations from the Orthogonality of 

the Gimbal Axis and Constraint Moment at Each Joint 

In this section, r-3 additional scalar equations are obtained by using the 
facts -that : 

(1) A. vector constraint torque at a typical joint can be isolated by 
summing the vector dynamical equations over all bodies that lie to 
one side of joint j 

(2) The constraint torque at joint j is orthogonal to the gimbal axis 
at joint j (i. e. , the dot product of the gimbal axis and the 
constraint torque at joint j vanishes). 

That is, for each gimbal axis i, if the equation 
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is summed over all bodies that lie to one side of joint i (that is, beyond gimbal 
axis g,, relative to body 0), then 


X £ i\ X V ' “h. = X'^ ^ + X e i>- X M \ 
\ \± \ J 


where 


€ =1 for all those bodies that lie to one side of joint i and is 0, 

i X 

otherwise. 

Note that the term 


X X < 

x JS J X J 



since the constraint torques cancel in pairs except at the joint j (when the 
constraint torques at joint j are summed over all bodies lying to one side of 
the joint relative to body 0). 

Using the fact that the constraint torque at the joint j is orthogonal to the 
gimbal axis at that joint, it follows that 


g; * 


^ i 6 ik ( ^ \ W |J. E > 

X \ (JL 


= 0 i = 1 , 2, . . . r-3 


(B- 10) 


Manipulation of Eq. (B-10) yields 


h • X ‘ix X k 


r-3 


+ 22 €v " gv 


|a X 0 / i 'k[Ji B k *k 

(JL \ k=l 


r-3 


(B- 11) 


i ^ 1 € ik ( E k 2 \ ^ \ 


£ k(J. g k \ = g 


k= 1 


.-x*. 


ik E k 
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In compact form, Eq. (B-ll) can be written as 


a. 


iO 


w 0 + 


r - d 

E 

Lr- 1 


a ik ^k 



whe re 


a i0 



6 i\ ® Xfj. 


a 


ik 



€ iX % ®k 


In Eq, 


(B- 12), 


a M is 
ik 


a scalar and a_ is 
iO 


a vector. 


(B- 12) 


B - 3 Set of r Dynamical Equations (Free From Constraint Torques) 

In this section, the set of r dynamical equations obtained in Section B-l 
and B-2 is given. Collecting the results of Eq. (B-7) and (B-12), it follows 
that 


r -3 

a oo ’ w o + E & 0k ^k = ^ E^ , 3 scalar equations 
k= 1 k 


r-3 


a. ^ I - 

iO 0 


X a ik^k = g i ' X £ ik E k’ r_3 e( i uations with 


(B- 13) 


k=l 


i = 1, 2, ... r-3 
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where 


‘00 = E XI % ' a dyadi< 

X [si 


‘Ok = E E V V ' g k- a vector 

X |JL 


E k = E x 


E 


r-3 


r-3 


'Xfx ^ v e k|ju g k ^k E X 


M- 


k=l 


EE 

k=l \ [x 


£ kjJL \(JL * g k / Y k 


a. 


EE 

X \i 


io = g i • 2 ^ 2 -j € a %’ a vector 


a ik = % • E E 

X jsX 


e iX e kn % • «k • a scalar 


e identifies those bodies u which sense the relative rate v. g- 

k(Ji ^ 1 k to k 

identifies the bodies X over which the dynamical equations are 
summed to isolate the constraint torque at joint i 

In matrix notation, Eq. (B-13) can be written as 


a 00 a 01 a 02 **• a 0, r-3 


a 1 0 a ll 


l,r-3 


a r-3,0 a r-3,l 


r-3, r-3 


/“ 0 


I ? 

*1 • E 


^ r -3 


E? 


£ 1 X E i 


g 


r-3 


E-,, 


E ' 
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It is understood that the dyadics and vectors appearing in Eqs. (B-7) 

and (B-12) are replaced by matrices of their components in Eq. (B-14). It 

is also noted from the definitions of a nrk , a A , , a_, a., that the matrix 

00 Ok lO lk 

involving these elements in Eq. (B-14) is symmetric. The dimensions of the 

matrix elements are as follows: a „ „ --3x3, a., --3x1, a, A -- 1 x3, 

00 Ok kO 

and a^ is a scalar for i, k = 1, 2, ... r-3. 

Equation (B-14) represents the desired set of r dynamical equations 
governing the motion of an n-hinged rigid-body spacecraft. 


52 


JPL Technical Memorandum 33-613 



JPL Technical Memorandum 33-613 


Table B-l. Set of r Dynamical Equations for an n-Hinged Rigid-Body Spacecraft 

(Free from Constraint Torques) 


Item 

Equations 

Remarks 

Starting point for 
obtaining set of 
n dynamical 
equations (free 
from constraints) 

Evv = 

P 

where 

E x = M x + D x x 

e x + E <■ 

j‘ J x J 

f x + 'y d xii x ^ o °xh 

P A p A 

Set of equations for n-hinged rigid-body spacecraft 
(involving constraint moments) that were derived in 
Appendix A; is constraint moment at joint j on 

body X; all other terms are defined in Table A-l. 


x 1 m u 

L H 

% x V 1 ! ■ “x x \\ ' w x + E M x° 

i‘ J X J 



\x = °x + x x 

- D x + ["\( D x u - D X D x) 



+ E "V 

p A 

,« u - W)] 



V = - m K 

X- % u - W] 


Elimination of 

constraint 

moments 

EE V" 

\ p 

since 

v - E E x 
\ 

Summing the constraint moments on body \ over all 
joints jcJ^ and over all bodies removes the constraint 
torques since they cancel in pairs (Newton s third 
law) 


EE 

x je j x J 

= 0 


Simplification 
of result 
obtained by 
summing 
over all 
bodies 

E Ev 

X p 

r-3 r-3 

u 0 + ^ > *kp + ^ 1 £ kp ^k ^k 

k=l k= 1 

g^ is the unit vector representing the gimbal axis at 

joint k; specifies the bodies p which sense the 

rotation v, g. ; is 1 if the gimbal axis is between 

’ k B k kp 

body p and body 0, otherwise = 




Table B-l. Set of r Dynamical Equations for an n-Hinged Rigid- Body Spacecraft 

(Free from Constraint Torques) (contd) 


Equations 


Elimination of 
constraint 
moments (contd. ) 


Simplification 
of result 
obtained by 
summarizing 
over all 
bodies X. 


= “o + E 


Rearranging, 


EEvvElEE 

\ p k= 1 ' X p 


p. € kp. ^k I ^k 


EvEvE 


f_j -Xti L-i V ' g k Y k ' 

X p k=l X 


E< 


Compact 

form 


L - J 

w 0 + Xj a ° k ^ k = Xj Ex 


aQQ is a dyadic; 

E^, E^ are vectors; result represents 3 scalar equations 


r-j 

Ev'§ v ■ 


.oTE v 

\ p 


Determination of 
additional r-3 
scalar equations 

Isolating con- 
straint moment 
at joint j 


vTE 

\ p 


^Xp kp ®k 


€ i\ 1 ^ , ^Xp «0 + ^ 1 e kp ^k ®k 
X p L k=l 




Summing the vector dynamical equations over all bodies that 
lie to one side of joint i isolates the constraint torque ; 
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Table B-l. Set of r Dynamical Equations for an n-Hinged Rigid- Body Spacecraft 

(Free from Constraint Torques) (contd) 
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Table B-l. Set of r Dynamical Equations for an n-Hinged Rigid-Body Spacecraft 

(Free from Constraint Torques) (contd) 


Item 

Equations 

Results 

Set of r dynamical 
equations for 
n -hinged rig id - 
body spacecraft 
(fre'e from con- 
straint moments) 

r i i / v f* \ 

a 00 a 01 “* a 0, r-3 w 0 Z-rf X. 

a 10 a ll “* a l,r-3 Vi g l 

a r -3, 0 a r~3, r-3 ^r-3| ^ <r - 3 \ 

Matrix components of dyadics and vectors are used 

in this equation; a^Q is 3x3, a^ is 3 x 1, a . ^ is 

1 x 3, a., is a scalar for i, k = 1, 2, ... r-3; 
lk 

matrix involving Sqq, a.^, a.^ is symmetric. 



APPENDIX C 


DETERMINATION OF EXPLICIT FORM FOR DYNAMICAL EQUATIONS 
FOR A SPECIFIC TOPOLOGICAL ARRANGEMENT OF 
5-HINGED RIGID BODIES 

In this section, the general results described in Appendix B are applied 
to a specific topological arrangement of 5 hinged bodies. In particular, results 
are obtained in this section for the arrangement of bodies shown in Fig. 5-1. 
Hence, the explicit form of the dynamical equations for the case in which 
n = 5 or less is obtained. 


C - 1 3 Scalar Equations Obtained by Summing Dynamical Equations Over All 

5 Bodies ~ 

Three scalar equations are obtained by summing the dynamical results 
for each body \ over all 5 bodies. This result is written as 



Since the constraint torques cancel in pairs, the term 


E 


E 



o 


(C-1) 


(C-2) 


It can be seen from Fig. 5-1 that the angular velocities of the bodies are 


% = “0 

W 1 = w 0 + h g i 

w 2 = w 0 + h g l + ^2 g 2 (C ‘ 3) 

u>3 = + Y3 S3 

“4 = “0 + h g 3 + h g 4 
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or, in general. 


CO = co A + 
fa 0 


% Y k g k 


where identifies which angular rates y^. g^_ are involved in the angular 

velocity of body fa . Effectively, e ^ = 1 if the gimbal axis g^ lies between 

body fa and body 0 and is 0 otherwise. Using Eq. (C-4) in Eq. (C-l), it 


follows that 


4 4 


4 4 


£ £ v •*.•£ £ v £ 


V g k^k 


X = 0 |a=0 


X=0 p =0 


EK‘ E v E % s k v k 

X= 0 \ p =0 k= 1 


Equation (C-5) can be rewritten as 


4 4 


4/4 4 


ZEvvEEE 


X=0 p =0 


k= 1 \ X = 0 p = 0 


£ kp % / g k ^k 


4/4 4 


EEE 

k=l \ X= 0 jjl =0 


£ kp % ) ‘ g k^k 


In compact form, Eq. (C-6) becomes 


a 00 • “0 + 2-J (b 0k • S k > V k = E x - Ej b 0k ' 8k i k 
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where 


4 4 

‘oo= EE 

\ = 0 p = 0 


% = b 00 


4 4 


} 0k ^ j ^ 1 € k[Ji ^p 

\ = 0 p =0 


a 0k ~ b 0k ’ g k 


The dyadics b^ for k = 1, 2, 3, 4 can be evaluated by inspection of Fig. 5-1 if 
the meaning of is kept in mind (recall that identifies those bodies p 


which sense the angular rate y^ g^_ ). 


It follows that 


4 2 

E E 

\ = 0 [i = l 


(since only bodies 1 and 2 
sense y x gj) 


4 4 


\-0 p =3 


(since only body 2 senses y^ g^) 


E E (since bodies 3 and 4 sense y^ g^) 


Y, \4 

\ = 0 


(since only body 4 senses y^ g ) 
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C-Z Equations Resulting from Orthogonality of Constraint Moments 
and Gimbal Axes 


In this section, the additional four equations resulting from the 
orthogonality of the constraint moment and the gimbal axis for each of the 
four joints are obtained. First, it is seen from Fig. 4-1 that the constraint 
moment at joint 1 is isolated by summing over bodies X. = 1, Z; that the con- 
straint moment at joint Z is isolated by considering body Z only; that the con- 
straint moment at joint 3 is isolated by summing over bodies 3, 4; and that 
the constraint moment at joint 4 is isolated by considering body 4 only. From 
the orthogonality of the constraint moment and the gimbal axis at each joint, 
it follows that 


g • = 0 

1 

g 2 * = 0 

g 3 * = 0 



(C-9) 


where a typical constraint moment M = M for X = 1, j = 1. 

j A 1 

For a fixed i, it follows that 


M x. = I] e ix X) M x. 
3 K s ej x 3 


with j = i 


«i * M x. = 0 

J 


(C- 1 0) 


where € identifies those bodies X over which the sum is taken. It has already- 
been observed that for i = 1 the sum is taken over X = 1, 2; for i = 2, only 
X = 2 is involved; for i = 3, the sum is taken over X = 3, 4; and for i = 4, only 
X = 4 is involved. Using Eq. (C-10), it follows that 
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(C-ll) 


X x. _ 2— 

J X 


® • u> - E 

Xp p X 


h ■ K. - 0 = «i • E V E % ' "n ' E X 

j X \ p=0 


Expansion of Eq. (C-ll) yields 


^ 1 \p w 0 + ^ j $ Xp ^ ; *kp g k ^k 

p =0 p =0 k=l 


Yu ^ 'kp ^ k ~ Ex 


Simplification of Eq. (C-12) yields 


vE 


e ix Y v ■ vsi- Y\Y* ix Y 


-i i x L-U 

X p =0 


k=l \ X 


h Y € ix E x- g i- Y 


k= 1 \ X 


k|Ji \\i 


/ ®k ^k 


In compact form, Eq. (C-13) can be written as 


h ' b i0 • S + g i ' Y (b ik ’ g k> % 


i • E 'iK E x - g i ' E (b ik ■ g k> il 


= % • 2 _! 'a E x - g i ' 2-i ( 
X k= 1 
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where 




b 


ik 


4 



X |JL =0 


€ kp 



In the notation of Appendix B, Eq. (C-14) can be written as 


a. 


iO 


CO *f 

o 


E 

k= 1 


a ik ^k 


= S t 


E e a E K-E 


(Si • h ; 


ik 


V \ 


(C- 15) 


k=l 


where 


a 


a 


iO 

g i • 

b i0 

ik = 

g i * 

b ik 


(C - 16) 


The dyadics can be written as (see Table C-l) 



2 4 


E E V 

X= 1 p = 0 


(since X. = 1, 2 are involved in sum over X.) 


b 


20 


4 


E 

p =0 



(since only X. = 2 is involved in sum over X) 


b 


30 


EE 


\=3 \x -0 



(since X = 3, 4 are involved in sum over \) 
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(since only X = 4 is involved in sum over X) 

(C- 17) 


b 


40 



Note that the bodies p, involved in the sum 


^ j € kfx \fji 


fJL =0 


for k = 1, 2, 3, 4 have already been identified (€ identifies the bodies p 
which sense the rotation g^). Now the bodies k which are picked out by 
the sum 


E 


€ . 


ik 


have also been identified for each i. Consequently, the bodies p, k involved 
in the sum 


b 


ik 





can be written down by inspection. Written in terms of an array, elements 

b™, b^, (derived in Section C-l), and b t b., become 
00' Ok kO ik 
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4 4 

E E % 

\=0 p .=0 

2 4 

EE % 

\=1 |JL =0 

E % 

p. =o 

4 4 

E E v 

V=3 p, =0 
4 

E v 

(JL =0 


EE v 

\ = 0 p. = l 

EE v 

\=1 p. — 1 

2 

52 $2 M- 
H = 1 

E E % 

\ = 3 (JL =1 

E V 

(L = l 




00 

b oi 

b 02 

b 03 

10 

b ll 

b 12 

b l3 

20 

b 21 

b 22 

b 23 

30 

b 31 

b 32 

b 33 

40 

b 41 

b 42 

b 43 


4 4 

E E v 

X = 0 (jl =3 

EE v 

\= 1 (JL =3 
4 

52 

M- =3 
4 4 



\=3 p. =3 
4 


E % 

[i-=3 




(C- 18) 


Note that Eq. (C-18) could have been written immediately by inspection of 
Fig. 5-1! In row 1 of the array, the summation over X is from X = 0 to X = 4; 
in row 2 (corresponding to the dot product of g^ and the constraint moment at 
joint 1) the summation over X is from X = 1 to 2; in row 3 (corresponding 
to the dot product of g^ and the constraint moment at joint 2), only X = 2 is 
involved; in row 4 (corresponding to the dot product of g^ and the constraint 
moment at joint 3), the summation over X is from X - 3 to X = 4; in row 5 
(corresponding to the dot product of g^ and the constraint vector at joint 4), 
only X = 4 is involved. Similarly, in column 1, the summation over jx is from 
0 to 4 (all bodies sense coq); in column 2, the summation over \i is from 
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1 to 2 (only bodies 1 and 2 sense y^ g^); in column 3, the summation over \x 
involves only p ~ 2 (only body 2 senses y^ g^ ““ ^ end of a chain! ); in 

column 4, the summation over \x is from 3 to 4 (only bodies 3 and 4 sense 
y 3 g 3 ); in column 5, the summation over p involves only body 4 -- it, too, is 
the end of a chain! ). 


Set of r (r = 7) Dynamical Equations for a 5-Hinged Rigid-Bod^ 
Spacecraft 


In the section, the set of 7 dynamical equations governing the behavior of 
the 5-hinged rigid-body spacecraft shown in Fig. 5-1 are given. These equa- 
tions are obtained by collecting the results generated in Sections C-l and 

C-2. 


In the following equation, it is understood that the matrix components of 
operations involving dyadics and vectors are used (recall, that the dyadic s b^^ 
with ft, m = 0, 1, . . 4 are defined in Eq. (C- 18)) 
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It is worth repeating the fact that the dimensions of the elements are as 
follows: 

b 00 - 3 x 3, b ok • g k - 3 x 1, g k • b k0 -- 1 x 3, g. • b. k . g k -- 1 X 1 
for i, k = 1 , ... 4. 

As mentioned previously, the square matrix on the LHS of Eq. (3-19) is sym- 
metric so only 


s = £ (a + St) = 4 (5 + 1) = 15 


of the 25 elements must be evaluated (note, however, that b^Q is a symmetric 
matrix, too!). Typical elements of this array are the 3x3 matrix b^, the 


3x1 vector b^ • g^ and the scalar g^ • b^ • g^. 
of this array were previously, defined as (see Eq, 

a oo ■ b oo 



3 Ok = b 0k • *k’ 


k = 1, 2, ... 4 

o 

II 

OQ 

• 

cr 

s* 4 

o 


k = 1 , 2, . . . 4 

a.. = g. • b., • 

lk l lk 

*k' 

i, k = 1 , 2, . . . 


(C-20) 


In compact form, Eq, (C-19) can be written as 

A co = L (C-21 ) 

where A is a 7 x 7 matrix, go is a 7 x 1 vector (consisting of oo^ and the relative 
angular accelerations y^), and L is a 7 x 1 vector (which can be considering the 
forcing function of the matrix differential equations). 
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C-3. 1 Transformation of Elements into Appropriate Coordinate System 


In this section, the vector basis in which the elements must be expressed 
in order for the operations to be performed is selected. A logical basis is that 
associated with body 0 (the base body). For the attitude control problem, it is 
usually necessary to control only the attitude of the base body (the attitude being 
kinematically related to the angular velocity vector cOq). Hence, the transfor- 
mations from bodies 1, 2, 3, 4 to body 0 are required. 

C-3. 1.1 Coordinate Transformations 

The transformations required in performing the operations depicted in 
Eq. (C-19) are (note these are not restricted to small angles) 


1 T - 

c 0 = cos Yj E + (1 - cos Yj) gj g i - S Yj gj 

C i = cos Y 2 E + (1 - cos v 2 ) g 2 §2 - S Y 2 g 2 

(C-22) 

3 T 

c 0 = cos Y 3 E + (1 - cos v 3 ) g 3 g 3 - S y 3 g 3 

C 3 = cos Y 4 E + (1 - cos y 3 ) g 4 g 1 - S y 4 g 4 

where 

E is the 3x3 identity matrix 
T 

g^ g^_ is an outer product, k = 1, 2, ... 4 

g^ is a 3 x 3 skew symmetric matrix representing the cross product 
operator x, k = 1, 2, ... 4 

k 

C k j ^- s transformation from the coordinate frame of body k-1 to 
that for body k 
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The results given in Eq. (C-22) were obtained in this analysis by 
manipulation of the relationship between two coordinate systems in terms of 
Euler’s parameters (for a different approach see Ref. [5]). Consequently, 
the transformations relating bodies 1, 2, 3, 4 to body 0 can be expressed as 
(see Fig. 4-1) 

C* = cl 

0 0 


cr = 


c z c 1 

c i c o 


(C-23) 


4 4 3 

C 0 = C 3 C 0 


Note that in Eq. (C-22), it is assumed that the components of g^ are known in the 
coordinates of body 0, the components of g^ are known in the coordinates of 
body 1, the components of g^ are known in the coordinates of body 2; and the 
components of g^ are known in the coordinates of body 3. 

Mbreover, the direction cosine matrix relating the coordinate system for 
body 0 and the coordinate system of the interial frame N for a 3-2-1 sequence of 
rotations is 



C 1 C 2 


S 1 C 2 

" S 2 

c° - 

C N " 

- S 1 C 3 + c 

1 S 2 S 3 

C 1 C 3 + S 1 S 2 S 3 

C 2 S 3 


C 1 S 2 C 3 + 

S 1 S 3 

S 1 S 2 C 3 - C 1 S 3 

C 2 C 3 

where S. = sin 0., C. = 
1 11 

cos 0. 

l 

(with i = 1, 2, 3) 



(C-24) 
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The attitude angles 0 , 0^, 0^ of body 0 are related to the angular 
velocity coq according to 


\ 


S 1 S Z 


C 1 S 2 


-S, 


00 . 


00 . 


\ 3 / 


(C-25) 


C-3.2 Evaluation of Elements b * and a^ 

— — Jim Jem 


The following conventions are used in evaluating the dyadics b^ m and a^ 
m = 0, 1, ... 4: 

(1) A typical element $ (X, p = 0, 1, . . . 4; X ^p) is evaluated by 

expressing all vectors involved in the coordinate system of body X 
(i. e. , the basis of the first subscript of ^ ); e. g. , 


m 


- m [D 12 • D 21 - D 12 D 21 ] 


where is computed in the coordinates of body 1 (the first 

subscript of D ) and then transformed to the coordinates of body 2 

\p 


by 


* D 12*2 " C 1 1 2b 


and * s a l rea( iy known in the coordinates of body 2; consequently 

[®2lh = - m [( C l D l 2 )' D Z1 - ( C l D 12) D 2ll 

(2) A typical element b^* ^ = 1, . . . 4 is evaluated by transforming 

all ’ s from the coordinate system to body X to that of body 0 
according to 

[ %Jo = C X [ %A c 0 ’ withX ^’ = 0, 1, ... 4, 
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e • § • f 


b 01 ^01 + $ 02^0 + C l ^ $ 1 1 + $ 12 ^ 1 C 0 + C 2 ^21 + ^ 22-^2 C 0 


+ C 3 ^31 + $ 32 ^ 3 C 0 + C 4 ^*41 + $ 42 ^ 2 C 0 


(3) A typical element (i, k = 1, 2, ... 4) is evaluated by transforming 
all 4^ *s involved to the coordinate system of body i, e.g. 

[b 13 ]i = [* 13 + ® J4 ] 1+ cJ [® 23 + * 24 J 2 C 2 


^33] 3 “ + + [ ®i** + ^ 


33 ' 34 J 3 1 4 L 43 ' ^44 J 4 ~3 


Similarly, the operations b^ • g^_ and g^ • b^ • g^ involved in the a ^ m ,s 
are evaluated according to 


(1) b 0k g k f b 0k^0 ' ( C k g k) 

(2) V b ik' g k = (g i>i • t b ik J i ‘ ( C k g k) 


These results are summarized in Tables C-l and C-2. 


C-3.3 Evaluation of Forcing Function L 

In this section, the terms involved in the evaluation of the 7x1 forcing 
function L (the RHS of Eq. C-19) are given (see Table C-3). The column 
matrix L can be partioned as 



where is a 3 x 1 vector and is associated with body 0 and is a 4 x 1 
vector associated with bodies 1 thru 4. 
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SD 

For notational convenience, the terms are redefined as t . and 

represent the spring -damper moment acting on body X at joint j. 


C-3.3,1 Viscoelastic Moments and Joints 

In this section, the assumed form of the spring-damper interaction 
moments acting on body X at joint j is given. Equations representing the 
interaction moments are given by 


T 12 " t 22 

T n = - K i Y i - B i h 

(C-26) 

T 01 = " T 11 

T 44 = ' K 4 Y 4 ' B 4 Y 4 

T 34 = " t 44 

t 33 = - K 3 ^3 - B 3 Y 3 

t 03 = ' T 33 


Equation (C-Zfe) implies that 
acting on bodies 0 thru 4 are 


the total interaction spring -damper torques 


T 01 + T 03 K 1 + B 1 h + K 3 V 3 + B 3 b 
T 11 + T 12 K 1 Y 1 ' B 1 Y 1 + K 2 Y 2 + B 2 Y 2 
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t 33 + T 34 “ “ K 3 Y 3 " B 3 h + K 4 Y 4 + B 4 ^4 
T 44 = " K 4 ^4 - B 4 Y 4 


(C-27) 


C-4 Explicit Equation for (The Angular Velocity of the Base Body 0) 

In this section, the explicit form for the equation governing a> 0 is given. 
This can be obtained by partitioning the equation 


A w = L 


and eliminating the relative angular velocities from the equation. In partitioned 
form, Eq. (C-21) becomes 


r- 

A 1 , 

1 A . _ 

1 1 

l 12 

A „ 

1 A„ 

21 

i 22 J 




(C-28) 


where the dimensions of A^, A^, A^ 
3x4, 4x3, 4x4, 3x1, 4x1, 3x1, 
Eq. (C-28) yields 


’ A 22' 
4x1, 


W 0’ L 0’ 
respectively. 


and L are 3 x 3 , 

XV 

Manipulation of 


A 11 “0 + A 12“r 

= L 0 

(C-29a) 

A 21 w 0 + A 22 w R 

= l r 

(C -29b) 


Solving Eq. (C-29b) for co- and substituting the result into Eq. (C-29a), it 

XV 

follows that 




[An “ A 12 A 22 A 21^ L 0 _ ^ A ll “ A 12 A 22 A 21^ A 12 A 22 L R 


W R A 22 L R “ A 22 A 21 “o 


(C-30) 
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Equation (C-30) suggests that the RHS of the equation for could first be 

evaluated; next, the expresion for can be used in evaluating the RHS of the 

. * * • 

equation for co.^; then, and co^ can be obtaining by integrating co^ and co^ , 

• • 

respectively. The fact that can be expressed in terms of co^ is significant 
when one is interested in arriving at a suboptimal Kalman filter in which the 
dimension of the state to be estimated is kept as low as possible. 


Table C-l. Expression for Dyadics f , m = 0, 1, ... 4 
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Table C-l. Expression for Dyadic s 6., m = 0, 1, ... 4 (contd) 


Item 

Equations 

b 03 

4 4 

YLYu V = [ ®03 + ®04b + C ? f*13 + *14ll C> 0 
\ = 0 p =3 

+ C 2 f®23 + ®24h C 0 + C 3 f®33 + *346 C 0 
+ C ! [* 43 + ®44l4 C 0 

b 04 

4 

= ^*04^0 + C 1 ^14^1 C 0 + C 2 ^24^2 C 0 

\ = 0 

+C 3 f $ 34^3 C 0 + C 4 f^44^4 C 0 

b ll 

2 2 

2 -J = f *11 + *12 1 1 + C 2 [®21 + $ 22 ]2 C 1 

X= 1 p=l 

b 12 

2 

^2 = [® 12 h +C 2 [* 22^2 C 1 

\=1 

b l3 

2 4 

Yj I] V * [*13 + ®14h +C 2 [®23 + ®24h C 1 

[1=1 |JL =3 

b 14 

2 

Y1 *U = t *14 1 1 +c 2 [® 2 4 h C 1 
X= 1 
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Table C-l. Expression for Dyadics l 7 , m = 0, 1, . . . 


4 (contd) 


Item 

Equations 

b 22 

[4- ] 

1 ^22 J 2 

b 23 

4 

E ®2|i = ^23 + ^24 ^ 2 
p =3 

b 24 

f $ 24^ 2 

b 33 

4 4 

EE [ ®33 + *54 >3 +c 4 [®43 + * 44)4 C t 

\-3 [x =3 

b 34 

4 

$ X4 = ^ *34 J 3 + C 4 f *44 ^ 4 C 3 

X=3 

b 44 

[ ^ 44 J 4 

b im 

(b m |) T , i, m = 0, 1, ... 4; Hi m 
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Table C-2. Expressions for Elements a^^ (H, m = 0, 1, • . • 4) 


Item 

Equations 

a oo = b oo 

f $ 00 + %1 + %2 + + %4 ] 0 


+ C 1 f®10 + ®ll + *12 + »13 + »14>l c l 


+ C 2 [®20 + ®21 + ®22 + ®23 + ®24 ] 2 C 0 


+ C 3 [*30 + $ 31 + $ 32 + %3 + %4 I 3 C 0 


+ C 4 [< f 40 + *41 + *42 + t '43 + *44l4 C 0 

a oi 

[ t $ 01 + ®02 ^0 + C 1 1 + ®12 ^ 1 C 0 + C 2 ^ ®21 + ®22 ^2 C 0 

+ c° [« 3 j + * 32 ] 3 c^ + c» [® 41 + <f 42 ] 4 c o] ' c i 

a 02 

[f $ 02^0 + C l ^®12 ^ 1 C 0 C 2 ^22^2 C 0 + C 3 ^®32^3 C 0 


+ C 4 t*42)4 C o] • C 2 <S 2 >2 

a 03 

[[*03 + *04]o + C ? t*13 + *14]l C J +C 2 [ $ 23 + ®24l 2 C 0 


+ C 3 [®33 + *34l3 C 0 +C 4 (^43 + *44 U C 0 1* C ? <«3>3 

a 04 

[[®04lo + C ?[®14ll C J +C 2 [®24l2 C 0 + C § ^34^3 C 0 


+ C 4 f*44) 4 C o]- C S<*4>4 
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Table C-Z. Expressions for Elements (H , m = 0, 1, ... 4) (contd) 


Item 

Equations 

a n 

( g l)l [b *11 + *12b + C 2 ^ *21 + *22^ 2 C 1 ] ’ 

a 12 

( g l)l [ f *12 b + C 2 ^*22 b C 1 ] ' C 2 (g 2^ 

a 13 

( g l)l [ ^*13 + *14 ^ 1 + C 2 f *23 + *24 b C 1 ] * C 3 * g 3^3 

a 14 

( g lh [ f *14 h + C 2 ^ * 24 ^ 2 C 1 ] * C 4 (g 4 U 

a 22 

{g 2 } 2 • ^ *22 b (g 2 } 2 

a 23 

^ z\z ' t *23 + *24 b * C 3 ^ g 3^3 

a 24 

^ g 2^2 * f* 24 ,b * C 4 ^ g 4^4 

a 33 

( g 3 ) 3 ' £ [$33 + $ 34 ] 3 + c 4 [$43 + $44 ] 4 C 3 ] (§3)3 

a 34 

( g 3 ) 3 * [ [$3413 + C 4 [ $44 ] 4 C 3 ] - C 4 (§4)4 

a 44 

(£4)4 * [^44] 4 * (§4)4 

a ko 

(.*)* k = 1 » z > • • • 4 

a fm 

( a m ^) T , i, k = 1 , 2 , ... 4 , i 4 k 
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Table C-3. Evaluation of Forcing Function L 
(Appearing in A ( 1 ) = L) 



E 0 M o 


+ D.. x F_ + / D n x f F + mu x (m x D _ )] 

0 0 / j Op p L p p p pO J 

p A 


" w 0 X *00 * W 0 + K 1 Y 1 + B 1 Y 1 + K 3 Y 3 + B 3 Y 3 


Ej M 


1 +D 1 xF,+ / D n xC^TF + mm x (m x D , ) 1 

111 / j lp p L p P P pl 

p A 


[( C 0 “0 + V] gj) X [* u ]j • (c* + Vj gj] 


+ K 2 Y 2 + B 2 y 2 - Kj Yj - Bj Yj 


E 2 M 2 


+ D x F 0 + / D xC^[F + mm x (m x D ) 1 

2 2 / j 2jjl |jl l (jl jjl \x \xZ j 

p A 


[ C 0 w 0 + C 1 Y 1 g l + Y 2 g 2 ] X f $ 22^2 ' [ C 0 W 0 + C l Y l g l +Y 2 g 2 ] 


- K 2 Y 2 - B 2 Y 2 



M., + D- x F_ + / D_ x FF + mm x (m x D ,)] 

3 3 3 / ^ 3 p p L p p p p 3 J 

P^ X 

■ [°o“o + ^3 S3J x t *33)3 • [ c 0 “0 + ^3 g 3 ] 


+ K 4 ^4 + B 4 Y 4 - K 3 ^ - B 3 V 3 
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Table C-3. Evaluation of Forcing Function L 
(Appearing in A (I) = L) (contd) 

Item 

Equations 

E 4 

M. + D. x F. + / D, x C 4 [F + mw x (w x D ,)1 

4 4 4 / i 4 jjl. ta *- (jl p ' p p 4 J 

p A 

- [ C 0 “0 + C 3 ^3 g 3 + ^ g 4 ] * [*44^4 ' [ C o“o + C 3 ^3 g 3 
+ Y 4 8 4 ] -K 4 K 4 -B 4 ; 4 

*1 

C 0 “0 x 6 1 

g 2 

( C 0 “0 + °1 n *l) X g 2 

h 

c 0 ^0 X g 3 

h 

( C 0 w 0 + C 3 V 3 83 ) X g 4 

L o 

4 4 

Y A - Y (b 0 k • 8 k ) V k = E 0 + cj Ej + c® E 2 4 C° E 3 

\=0 k=l 

+ C 4 E 4 - ( b 01 ' C 1 81 ) Vj - (b 02 ' C 2 *z)Vi - ( b 03 ' C 3 83 ) V 3 
- ( b Q4 C 4 8 4 )v 4 
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Table C-3. Evaluation of Forcing Function L 
(Appearing in A (i = L) (contd) 


Item 

Equations 

L i 

g l ' ( E 1 + C 2 E z) 


" *1 ' f b l 1 ' *i ■»! + b 12 • ®2 ^2 + b 13 ' % V 3 + b 14 • g 4 ; 4 ] 

L 2 

g 2 • E z - S 2 ' t b 21 ’ Si V! + b 22 • g 2 y 2 + b 23 ' g 3 V 3 


+ b 24 ' ^4 J 

L 3 

g 3 ' ( E 3 + C 4 E 4^ " g 3 ’ ^ b 31 ’ g l h + b 32 ’ g 2 ^2 


+ b 33 * g 3 Y 3 + b 34 • g 4 V 4 ] 

L 4 

®4 ' E 4 - g 4 ' t b 41 ’ S; Y, + b 42 • g 2 Y 2 + b 43 • gj Yj 


+ b 44 ' g 4 Y 4 ] 
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